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Abstract 

A  macroscopic  electrochemical-mechanical  model  was  developed  to  predict  the  various  power  demands  of  an  electric  scooter  and  other  outputs, 
such  as  the  Li-ion  battery  pack  current  and  voltage  requirements  for  a  randomly  generated  drive  cycle  using  the  Matlab-based  Simulink  software. 
The  simulation  results  obtained  were  compared  with  the  actual  field  test  results  of  a  commercial  lead-acid  battery  pack  used  in  a  specific  electric 
scooter.  The  simulation  results  appeared  to  give  a  realistic  idea  of  the  dynamic  performance  of  a  conceptual  Li-ion  battery  pack  in  a  typical  drive 
pattern.  These  results,  though,  need  to  be  validated  in  an  actual  field  test. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Prior  modeling  efforts  in  the  Li-ion  battery  industry  have 
created  a  shift  in  focus  towards  developing  a  thermal  or  elec¬ 
trochemical  model,  or  a  combination  of  both  in  predicting  the 
thermal  and  electrochemical  response  of  the  Li-ion  battery  [1-4] . 
Mechanical  models  have  also  been  developed  to  calculate  the 
power  demands  of  the  battery  in  a  specific  drive  cycle  [5,6]. 
No  specific  works  have  been  reported  to  study  the  thermal  and 
electrochemical  performance  of  the  Li-ion  battery  in  a  scooter 
application  taking  into  consideration  the  drive  cycle,  mechanical 
parameters  of  the  vehicle  and  driving  conditions. 

This  modeling  effort  will  help  to  design  and  develop  Li-ion 
batteries  that  employ  a  phase  change  material  (PCM)  thermal 
management  system  for  an  electric  scooter.  Al-Hallaj  et  al.  [7] 
simulated  the  thermal  performance  of  a  Li-ion  battery  for  an 
electric  car  with  a  passive  thermal  management  system  using 
phase  change  material  to  demonstrate  its  feasibility.  In  an  elec¬ 
tric  scooter  application  where  space  constraint  limits  the  use 
of  active  thermal  management  system,  a  passive  PCM  system 
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proves  a  feasible  choice.  Khateeb  et  al.  [8]  described  a  detailed 
procedure  of  designing  Li-ion  batteries  for  an  electric  scooter, 
and  presented  simulation  results  showing  the  feasibility  of  a 
novel  passive  thermal  management  system  for  Li-ion  batteries 
using  phase  change  material  (solid-liquid).  Detailed  experimen¬ 
tal  results  were  presented  showing  the  thermal  response  of  the 
Li-ion  battery  in  the  presence  of  the  PCM  thermal  manage¬ 
ment  system.  Experimental  results  were  also  modeled  using 
a  simplified  lumped  thermal  model  [9].  It  was  important  to 
model  the  dynamic  electrochemical  and  thermal  characteris¬ 
tics  of  the  Li-ion  battery  during  actual  field  tests  as  a  function 
of  the  drive  cycle,  mechanical  parameters  of  the  vehicle  and 
other  parameters.  This  work  presents  the  development  of  an 
electrochemical-mechanical  model  to  evaluate  the  dynamic  cur¬ 
rent,  voltage,  state  of  charge  (SOC),  response  of  the  Li-ion 
battery  as  a  function  of  a  specific  drive  cycle  and  mechanical 
parameters  of  the  electric  scooter. 

2.  Model  development 

This  modeling  work  allows  for  simulation  of  the  performance 
of  the  Li-ion  battery  in  the  actual  riding  conditions  of  the  Zappy 
electric  scooter.  The  simulation  was  performed  using  Simulink 
(Matlab  Version  12.1).  The  mechanical  model  calculates  the 
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Nomenclature 

a  acceleration  of  electric  scooter  (ms-2) 

Ca  dimensionless  aerodynamic  drag  coefficient 

Ct  dimensionless  rolling  resistance  of  tire 

d  density  of  air  (kg  m-3) 

g  acceleration  of  gravity  (ms-2) 

Ip  combined  current  of  number  of  cells  or  string  of 
series  cells  connected  in  parallel  (A) 
mr  mass  of  rider  (kg) 

ms  mass  of  electric  scooter  (kg) 

rip  number  of  cells  or  string  of  series  cells  in  parallel 

ns  number  of  cells  or  string  of  series  cells  in  series 

Pace  power  demand  for  acceleration  of  electric  scooter 
(W) 

Aero  power  demand  to  overcome  aerodynamic  drag 
(W) 

Paux  power  demanded  by  auxiliary  systems — head 
lights,  signal  lights,  horn,  etc.  (W) 

^grade  power  demand  on  an  incline  (W) 

Ppar  parasitic  power  demanded  by  battery  system  (W) 
Ptotal  total  power  requirement  of  electric  motor  (W) 
^traction  power  demand  for  traction  of  vehicle  (W) 

Ptyre  power  consumed  for  motion  of  tires  (W) 

R  overall  resistance  of  Li-ion  battery  (£2) 

v  velocity  (speed)  of  electric  scooter  (ms-1) 

V  voltage  demand  of  Li-ion  battery  (V) 

Vo  open  circuit  voltage  of  Li-ion  battery  (V) 

Vs  combined  voltage  of  number  of  cells  in  series  (V) 

WQ\  electric  work  (W) 

Greek  letters 

r]m  efficiency  of  electric  motor  and  controller  system 
(75-80%) 

0  angle  of  inclination  of  road 


power  requirements  of  the  electric  scooter  as  a  function  of  the 
rider  weight,  acceleration,  road  inclination,  etc.,  whereas  the 
electrochemical  model  outputs  the  current  and  the  voltage  pro¬ 
files  from  the  Li-ion  battery  based  on  the  power  demand  from 
the  electric  motor.  The  voltage  profile  is  generated  in  the  simu¬ 
lation  via  an  integrated  empirical  plot  of  voltage  versus  state  of 
charge. 

Thus,  at  each  time  step  in  the  simulation,  the  capacity  of  the 
battery  available  is  calculated  based  on  the  current  consumption. 
As  a  next  step,  the  model  calculates  the  SOC  and  the  calculated 
value  is  compared  to  the  open  circuit  voltage  versus  SOC  plot. 
Then,  the  computed  battery  OCV  is  used  to  estimate  the  battery 
voltage  during  discharge.  The  battery  is  at  a  complete  charged 
state  when  SOC  is  100%  and  at  a  complete  discharged  state 
when  SOC  is  0%. 

Thus,  the  electrochemical  and  mechanical  model  is  integrated 
into  the  Simulink  and  simulation  is  conducted  under  various 
operating  conditions.  The  input  to  the  model  is  the  speed  versus 
time  data,  and  this  data  is  randomly  generated. 


The  important  features  of  this  electrochemical-mechanical 

model  are  as  follows: 

1 .  The  total  power  requirement  of  the  electric  motor,  the  power 
output  from  the  battery  and  the  state  of  charge  of  the  battery 
can  be  determined  as  a  function  of  the  speed  of  the  scooter, 
acceleration,  road  inclination  (grade)  and  rider  weight. 

2.  The  model  can  simulate  any  drive  cycle  required  and  verify  if 
the  Li-ion  battery  can  meet  the  power  requirements.  If  not,  the 
number  of  cells  needed  to  be  connected  in  series  or  parallel 
to  compensate  the  power  deficit  can  be  determined. 

3.  Safety  circuits  have  become  an  integral  part  of  Li- 
ion  batteries  for  high  wattage  requirements  of  any 
application — especially  automotive  requirements.  The  safety 
circuit  feature  can  also  be  implemented  in  this  model  where 
the  current  drawn  from  the  battery  will  be  limited  after  a 
threshold  limit  is  reached. 

2.7.  Assumptions 

1.  The  empirical  data  of  open  circuit  voltage  versus  SOC  used 
in  the  Simulink  to  predict  the  voltage  profiles  was  obtained 
from  C/1  rate  (1-h  discharge  rate)  constant  current  discharge 
test,  though  the  battery  could  discharge  at  varying  rates  in 
the  actual  drive  cycle. 

2.  The  resistance  versus  DOD  data  is  obtained  from  Sony  1 8650 
(1.8  Ah)  Li-ion  cell  at  C/1  rate  and  the  open  circuit  voltage 
(OCV)  data  is  obtained  from  Sony  1 8650  (1.4  Ah)  at  C/1  rate, 
whereas  Moli  energy  18650  (2.2  Ah)  Li-ion  cells  are  used  in 
the  experimental  work. 

2.2.  Traction  power  demand 

The  power  demand  for  traction  of  the  vehicle  can  be 


expressed  as  follows  [5,6]: 

Araction  —  ^grade  T  Acc  T  Ay  re  T  Aero  (1) 

Tgrade  =  ™gV  sin#  (2) 

Acc  =  mav  (3) 

Ayre  —  mgvC[  (4) 

Paero  =  0.5dCaAv 3  (5) 

Araction  =  mgv  sin#  +  mav  +  mgvC[  +  0.5dCaAv  (6) 

where  m  =  mY  +  ms. 

2.3.  Power  demand  of  the  electric  motor 


The  total  power  requirement  of  the  electric  motor  is  given  by: 
Aotal  =  +  Paux  +  Ppar  (7) 

7  in 

The  Li-ion  battery  is  designed  for  this  electric  scooter  using  a 
passive  thermal  management  system,  thus  eliminating  any  form 
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of  active  cooling  components  such  as  blowers  or  fans.  As  a  result, 
Ppar  =  0.  Also,  there  are  no  auxiliary  systems,  hence  Paux  =  0. 

2.4.  Power  supplied  by  the  Li-Ion  battery 

The  voltage  of  the  Li-ion  battery  is  calculated  using  the  fol¬ 


lowing  simplified  equation  [5,10]: 

V  =  V0 +  IR  (8) 

The  power  supplied  by  the  Li-ion  battery,  Pbatt  =  Vs/p.  (9) 

where 

Vs  =  nsV  (10) 

4  =  nvI  (11) 


The  discharge  as  a  function  of  depth  of  discharge  (DOD) 
was  obtained  for  Sony  (1.8  Ah)  18650  Li-ion  cell  from  litera¬ 
ture  [11].  The  open  circuit  voltage  data  was  obtained  for  Sony 
(1.4  Ah)  1 8650  Li-ion  cell  from  experiments  conducted  by  other 
members  in  our  research  group.  Both  of  these  empirical  data 
correspond  to  the  C/1  rate  discharge.  These  empirical  plots,  as 
shown  in  Figs.  1  and  2,  are  integrated  into  the  Simulink  model 
and  the  discharge  voltage  of  the  Li-ion  battery  can  be  generated 
as  a  function  of  depth  of  discharge. 

The  mechanical  and  electrochemical  model  is  implemented 
in  the  Simulink  as  shown  in  Fig.  3.  The  first  part  of  the  model 
block  diagram  represents  the  mechanical  model,  which  calcu¬ 
lates  the  total  power  demanded  by  the  electric  motor  from  accel¬ 
eration,  on  incline,  cruising  and  aerodynamic  power  demands. 
The  speed,  rider  and  vehicle  mass,  road  grade  and  other  scooter 
parameters  are  input  in  this  model. 

The  electrochemical  model  calculates  the  voltage  and  current 
output  of  the  Li-ion  battery.  The  input  to  this  model  consists  of 
the  empirical  results  of  the  open  circuit  voltage  and  resistance 

Resistance  of  Sony  18650  (1.8Ah)  Li-ion  cell  @  C/1  rate 


Fig.  1.  Discharge  resistance  vs.  depth  of  discharge  (DOD)  for  Sony  (1.8  Ah) 
18650  Li-ion  cell  [11]. 


Open  circuit  voltage  of  Sony  18650  (1.8  Ah)  Li-ion  cell  at  C/1  rate 


Depth  of  Discharge  (DOD) 

Fig.  2.  Open  circuit  voltage  (OCV)  vs.  depth  of  discharge  (DOD)  of  Sony 
(1.4  Ah)  18650  Li-ion  cell. 

as  a  function  of  depth  of  discharge.  The  sequence  of  simulation 
performed  in  the  Simulink  is  as  follows: 

•  Mechanical  model 

1 .  The  speed  data  from  a  driving  cycle  is  input  into  the  model 
along  with  other  parameters  specific  to  the  electric  scooter, 
rider  mass  and  inclination  of  road. 

2.  The  total  power  demanded  by  the  electric  motor  is  calcu¬ 
lated  as  a  summation  of  different  power  terms  according 
to  Eq.  (7). 

•  Electrochemical  model 

1.  Based  on  the  power  demand  from  the  electric  motor, 
the  current  requirement  of  the  Li-ion  battery  to  meet  the 
demand  is  calculated  according  to  Eq.  (9).  The  actual 
amount  of  current  that  the  Li-ion  battery  can  supply  is 
also  calculated  in  the  next  step  by  implementing  a  safety 
current  limiter  in  the  model  to  simulate  the  function  of  a 
safety  circuit. 

2.  The  electrical  capacity  consumed  of  the  Li-ion  battery  is 
calculated  by  time  integration  of  the  current  for  every  sim¬ 
ulated  time  step. 

3.  The  depth  of  discharge  (DOD)  of  the  Li-ion  battery  is 
calculated  by  dividing  electrical  capacity  consumed  by  the 
actual  initial  capacity  of  the  battery. 

4.  Based  on  this  DOD  at  every  time  step,  the  resistance  and 
open  circuit  voltage  of  the  battery  is  calculated  from  the 
empirical  results. 

5.  The  voltage  of  the  Li-ion  battery  on  load  is  then  calculated 
according  to  Eq.  (8). 

6.  The  above  steps  are  executed  in  the  loop  for  every  simula¬ 
tion  time  step,  which  generates  the  instantaneous  power, 
current  and  voltage  of  the  electric  motor  and  Li-ion  battery. 

3.  Simulation  results  and  discussion 

The  drive  cycle  shown  in  Fig.  4  was  randomly  generated  and 
used  in  the  dynamic  simulation  of  the  Li-ion  battery  module  in 
the  Zappy  electric  scooter.  The  various  inputs  to  the  simulation 
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MECHANICAL  MODEL 

Simulation  of  Li-ion  Battery  in  a  Zappy  Electric  Scooter 


Vo_disp 

ELECTROCHEMICAL  MODEL 


Fig.  3.  Simulink  block  diagram  of  mechanical-electrochemical  model  of  Li-ion  battery  in  Zappy  electric  scooter. 


Drive  cycle 
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Fig.  4.  Drive  cycle  simulated  in  Simulink  mechanical-electrochemical  model. 


according  to  Eqs.  (2-5)  are  summarized  in  Table  1.  The  drive 
cycle  data,  along  with  these  input  parameters  are  fed  into  the 
Simulink  model  and  the  simulation  was  performed  for  a  time 
period  of  4800  s.  The  results  from  the  simulation  yield  the  fol¬ 
lowing  important  information: 


Table  1 

List  of  values  of  parameters  used  in  simulation 


Parameters  Values 


Weight  of  scooter,  ras  (kg)  20 

Weight  of  rider,  rar  (kg)  60 

Density  of  air,  d  (kg  m-3)  1 .293 

Frontal  cross-section  area  of  the  electric  scooter,  Af  (m2)  0.5 

Acceleration  of  gravity,  g(ms-2)  9.81 

Dimensionless  rolling  resistance  of  tire,  Ct  0.008 

Dimensionless  aerodynamic  drag  coefficient,  Ca  1.2 


1 .  The  power  requirement  or  demand  from  the  electric  motor 
as  the  scooter  traces  through  the  drive  cycle. 

2.  Based  on  the  power  requirement,  the  discharge  current  output 
of  the  battery  can  be  determined  and  also  whether  the  battery 
will  be  able  to  sustain  the  current  requirement  throughout  the 
drive  cycle  period. 

3.  The  discharge  voltage  of  the  battery  and  the  time  duration 
for  the  battery  to  reach  the  cut-off  voltage 

4.  The  depth  of  discharge  (DOD)  of  the  Li-ion  battery. 

The  total  power  demand  from  the  electric  motor  is  shown 
in  Fig.  5.  The  electric  motor  is  rated  at  0.5  hp  (373  W)  power 
output,  though  the  motor  can  provide  short  bursts  of  peak  power 


Total  Power  requirement  of  the  Electric  motor 

500| — — — i  |  |  | — rrr —  j  i|  - n — 1 1  — ITT] - - 

450- 


Time  (s) 


Fig.  5.  Simulated  total  power  requirement  of  electric  motor  in  Zappy  electric 
scooter. 
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Current  output  from  the  Li-ion  battery  module 

40 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 


Time  (s) 

Fig.  6.  Simulated  Li-ion  battery  current  output  in  Zappy  electric  scooter. 

exceeding  the  power  rating.  The  average  power  requirement  of 
the  electric  motor  is  175-225  W. 

In  response  to  the  above  power  demand  from  the  electric 
motor  for  traction  of  the  electric  scooter,  the  discharge  current 
output  of  the  Li-ion  battery  is  shown  in  Fig.  6.  The  average 
current  throughout  the  entire  drive  cycle  period  is  about  15-1 8  A. 
The  average  current  output  of  the  lead-acid  battery  was  about 
1 8-24  A  on  real  time  field-testing  of  the  electric  scooter  with  a 
rider  weight  of  65  kg. 

The  actual  field  test  results  of  the  electric  scooter  using 
the  lead-acid  battery  is  shown  in  Fig.  7.  The  current  output 
of  the  lead-acid  battery  during  this  field  test  suggests  that  the 
battery  was  discharged  for  more  than  a  1-h  discharge  rate.  Ide¬ 
ally,  the  battery  should  be  discharged  within  2770  and  3600  s 
as  per  the  rated  capacity  of  18  Ah  but  during  the  field  test, 
the  18  Ah,  12  V  lead-acid  battery  was  able  to  discharge  for 
1950  s  before  the  electric  vehicle  halted.  Although,  the  com¬ 
mercial  lead-acid  battery  was  rated  at  1 8  Ah,  actual  laboratory 


Voltage  output  of  the  Li-ion  battery 

1 6 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - r 
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Battery  cut-off  voltage  reached 

80  500  1000  1500  2000  2500  3000  3500  4000  4500  5000 

Time  (s) 

Fig.  8.  Simulated  Li-ion  battery  voltage  output  in  Zappy  electric  scooter. 

tests  at  5  h  and  1-h  discharge  rate  delivered  less  than  the  rated 
capacity. 

In  the  case  of  the  simulation,  the  Li-ion  battery  with  a  capac¬ 
ity  of  24  Ah,  12.6  V  was  able  to  sustain  for  close  to  4400  s  before 
the  battery  is  not  able  to  provide  further  current.  The  average  cur¬ 
rent  from  the  Li-ion  battery  in  the  simulation  was  about  1 5-1 8  A, 
which  is  about  75%  of  a  1-h  discharge  rate  that  should  discharge 
the  battery  in  4800  s  as  per  the  rate  capacity.  One  of  the  reasons 
could  be  that  the  resistance  values  used  in  the  simulation  corre¬ 
sponded  to  a  1-h  discharge  rate  and  the  resistance  of  the  battery 
should  be  low  at  75%  of  the  1-h  discharge  rate. 

The  voltage  response  of  the  Li-ion  battery  is  shown  in  Fig.  8. 
The  average  voltage  of  the  Li-ion  battery  in  the  entire  drive 
cycle  is  found  to  be  12.5V  and  the  nominal  voltage  is  11.1  V, 
though  there  exists  some  peak  voltage  reaching  a  maximum  of 
14  V.  It  is  also  important  that  the  Li-ion  battery  does  not  over- 


Current,  Voltage,  Temperature  response  of  Lead-acid  Battery  in  Zappy 

Electric  Scooter 
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Fig.  7.  Preliminary  actual  results  of  Zappy  electric  scooter  with  a  commercial  lead-acid  battery  [8]. 
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Depth  of  Discharge  of  Li-ion  Battery  in  a  Zappy  Electric  Scooter 


Time  (s) 


Fig.  9.  Simulated  depth  of  discharge  (DOD)  of  Li-ion  battery  in  Zappy  electric 
scooter. 


charge  to  higher  voltage.  The  battery  is  found  to  reach  the  cut-off 
voltage  of  9  V  at  4400  s  and  the  battery  is  simulated  to  shut 
down  the  power.  The  depth  of  discharge  of  the  Li-ion  battery 
reaches  72%  at  the  cut-off  voltage  of  the  battery,  as  shown  in 
Fig.  9. 

4.  Conclusions 

The  work  presented  in  this  paper  described  the  development 
of  a  mechanical-electrochemical  model  of  Li-ion  battery  for 
an  electric  scooter  application.  The  Li-ion  battery  powering  a 
Zappy  electric  scooter  was  modeled  and  its  electrochemical  per¬ 
formance  simulated  during  a  randomly  generated  drive  cycle, 
by  means  of  a  specially  developed  electrochemical-mechanical 
model.  The  simulation  results  closely  predict  the  power  demands 
of  the  electric  motor,  instantaneous  Li-ion  battery  pack  voltage 
and  current  outputs. 

The  average  current  output  of  the  Li-ion  battery  in  the  drive 
cycle  was  about  15-1 8  A,  compared  to  1 8-24  A  measured  under 
actual  field  test  conditions  with  the  lead-acid  battery  powering 
the  scooter.  The  lower  current  demand  on  the  Li-ion  battery 
may  be  attributed  to  the  higher  average  voltage  output  of  the 
Li-ion  cells,  as  well  as  the  higher  discharge  capacity  compared 
to  the  lead-acid  battery.  The  discharge  capacity  of  the  Li-ion 
battery  was  24  Ah,  compared  to  1 8  Ah  for  the  lead-acid  battery. 


The  average  voltage  output  of  the  Li-ion  battery  was  around 
12.5V  (with  the  actual  nominal  voltage  of  the  Li-ion  battery 
being  14.4  V).  Once  the  output  voltage  would  hit  the  set  value 
of  9  V,  the  simulation  cuts  off  the  discharge  and  the  vehicle  would 
stop. 

Overall,  the  mechanical-electrochemical  model  appears  to 
give  a  realistic  prediction  of  dynamic  performance,  which  qual¬ 
itatively  agrees  well  with  expectations.  These  simulation  results 
have  to  be  validated  by  actual  field  tests  before  further  refine¬ 
ments  to  the  model  are  justified. 
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